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Abstract 

Model independent, mi, adiabatic, m_i, and high-energy, 771,3, energy weighted sums for the 
isoscalar and isovector nuclear excitations are investigated within the framework of the kinetic 
theory adopted to the description of a two-component nuclear Fermi-liquid. For both the adiabatic 
and scaling approaches, the connection of the EWS m_i and 7713 to the nuclear stiffness coefficients 
and the first- and zero-sound velocity is established. We study the enhancement factor ki in the 
energy weighted sum m! x for the isovector excitations and provide the reasonable explanation of the 
experimental exceeding of the 100% exhaustion of sum m\ for the isovector giant dipole resonances. 
We show the dependence of the enhancement factor kj on the nuclear mass number A and analyse 
its dependence on the Landau's isovector amplitude F[. 

PACS numbers: 21.60.Ev, 24.30.Cz 
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I. INTRODUCTION 



m k 



The strength function is the basic characteristic which determines the behavior of a 
quantum system in an external periodic field U ext = X(t)q (\{t) = X e~' lujt + c.c, here q is 
the transition operator) 

S(E) = ^ |(tt n | q \V )\ 2 6(E - E n ) 7 E = f)w, (1) 

where \l/ n and E n are the eigenf unctions and the eigenenergies of the total hamiltonian H, 
respectively. Using the strength function S(E) one can calculate the moments m k (EWS) 

I dES{E)E k = Y,\^n\ q \^ )\ 2 (E n -E ) k . (2) 

Here, for convenience, we have included the ground state energy E Q = into the energy 
factor. Special role of the EWS m k is caused by its connection to the transport characteristic 
of the system. For example, the sums m_i and ^-3 determine the stiffness and mass 
coefficients for the collective excitations in the system Determined via the properties of 
the ground state of the system the sum mi plays a specific role. In many cases, it does not 
depend on the model used for the description of the collective motion. This allows one to 
test the results of theoretical calculations as well as the correctness and the completeness of 
the experimental data. 

During a few years, significant attention was paid to the analysis of EWS for the giant 



multipole resonances 



, signmcant attention was paid to tne analysi 

(gmr) q, y, q, b, y, 0, y, y, [lo| . Nuclear gi 



iant resonances exhaust 



a significant part of EWS (sometimes near 100%) and establish the relatively simple con- 
nection between the values of and the basic characteristics of the GMR. However, some 
problems occur while researching the EWS for the isovector giant dipole resonances which 
are the best investigated experimentally. The connected problem is that the sum mi is not 
model independent because of the dependence of the effective nuclear forces on the nucleon 
velocity. Thus, for the theory to agree with the experimental data, one has to include a 
phenomenological enhancement factor to the sum mi pj. As a consequence, this leads to 
the modification of other sums m^ and can affect the definition of the nuclear transport 
characteristics. 

In this work, we study the EWS for the isovector collective excitations in heavy nuclei 
and nuclear matter. Our approach is based on the kinetic Landau- Vlasov's theory adopted 



to a two-component nuclear Fermi-liquid. In Section 2, we consider the connection between 
the EWS m,k and the linear response function. The connection of the Landau's theory of 
the Fermi-liquid to the hydrodynamical model and to the scaling approximation is shown in 
Section 3 [LlIQ]- m Section 4, we apply our approach to finite nuclei. The main conclusions 
of the work are formulated in Section 5. 



II. LINEAR RESPONSE FUNCTION AND EWS 

Let us consider the response of a nucleus on an external field U ext (t) periodic in time 
which is switched on adiabatically at t = — oo: 

U ext (t) = Xoe-^+^q + X^-^q*, (3) 

where q is the Hermitian operator, 

A 

q = ^2q(n,Tj), (4) 

i=l 

A is the mass number, and Tj is the isotopic variable. If Ao << 1, then quantum mechanical 
expectation of the operator q takes the following form (see. ["Kjl ]) 

(ff>=x(w)A e-*- t + x»ASe*- t , (5) 

where x{ u ) is the linear response function 

X^) = Y J \^n\q\^ )\ 2 
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E n — E — fojj — iO E n — E + + iO 

it' 

Let us introduce the polarization response function 



(6) 



It is easy to establish the connection between the EWS and the linear response 
function xi^)- Let us take the Taylor expansion of the function x^K^) i n a series in hw as 
u) — > (adiabatic expansion) and in a series in (fiw)' 1 as u — > oo (high-frequency expansion). 
Using ([2]) and ([7]), we have 

X W ML = 2 [m_i + (/^) 2 m_ 3 + ...] , (8) 
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[mi + {fiw) 2 "^3 + ••■] 



(9) 



Below we will pay a special attention to the investigation of the sums m_i, ni\ and 777.3. 
Using these sums, one can define two averaged energies of collective motion 



£1 



mi 
m_i 



E% 



7773 
777i 



(10) 



It is easy to see that the closeness of the energies E\ and £3 to each other determines the 
exhaustion of the EWS nik by one state \l/ n (see (T5])). If the effective nuclear forces do not 
depend on the nucleon velocity, then the sum nil can be easily calculated and takes the form 
which does not depend on the model of collective motion. Namely, 



777i 



q,H 



f ) = ^|(^|g|f )| 2 (K-i?o) 



2m 



— * (■ 

dfp eq (r)\Vq(r)Y 



(11) 



where p eg (r) is the nucleon density for the ground state of the nucleus 

t A 

/>,,,{>■) = { 



i=l 



r — Tj 



(Here, and in the following, the symbol "eq" means that the proper value is related to 
the equilibrium (basic) state of the nucleus.) The expression ffTT]) is the so-called model 
independent EWS rule. If only one (collective) state ^ n =G exhausts the sum rule (ELD), i.e., 



777i 



\(^ G \ q \^o)\ 2 (E G -E ), 



(12) 



then we have E x w E 3 from ([2]), (fIU|) . 

The low-frequency (adiabatic) sum m_i is connected to the nuclear stiffness under the 
adiabatic slow deformation of the nucleus, in another words under the deformation that does 
not lead to the quantum transitions between nuclear levels. To reveal this connection, we 
will evaluate the energy variation AE of the nuclear ground state in an external static field 
Uext = Ao<? for Ao — > 0. Using the quantum perturbation theory for the calculation of the 
wave function \1/ of Hamiltonian H' = H + Ao<? in the second order in the small parameter 
Ao, we obtain 

AE ad = (tf I H - (*o| H |^o) = Agm_ x . (13) 
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Let us calculate the variation of the nuclear form parameter Q = (^\q\^) in the external 
field Xoq, 

AQ = Q = (tf | q |tf ) - (tf o| q l^o) = 2 A m_ 1; (14) 

where we have assumed (^q| q \^q) = 0. From (1131) and (I14p . we find the nuclear stiffness 
parameter Cq A cL with respect to the adiabatic change of the nuclear form as 

d 2 AE ad 1 

C Q , ad - - — . (15) 

Let us now consider the high-frequency sum m%. We introduce the wave function ty sc , 
which is obtained from the wave function of the nuclear ground state \l/o by means of the 
scale transformation (scaling-approach), 

# sc = e"fr«* , (16) 

where v is the small parameter of the scale transformation. 

In the case of a many-particle wave function given by the determinant built on the 
one-particle wave functions Q (r), the exponential operator of the scale transformation in 
( |T6l) acts on each function 4> a (r) independently For example, at the quadrupole deformation 

$ = £(r?-3z?), 
i=i 

one can see from (Tl6|) that ^ sc is also a determinant which is built on the one-particle 
functions (p a)SC (r) obtained by the scale transformation of coordinates. Namely, 

(j> a ,sc(f) = (/) a>sc (x,y,z) = (j) a (e u x, e u y, e~ 2u z), (17) 

where v = —2h 2 v/m. As can be seen from (ITTj) . the scale transformation does not violate 
the orthonormalization of the wave functions. By means of (fTEl) . the energy change AE can 
be found within the scaling approximation as 

AE = (* 8C | H |* sc ) - (* | H |^ ) = v 2 m 3 . (18) 

Using ( fl6|) . we obtain the connection between the parameter of scale transformation v and 
the deformation parameter Q: 

Q = (* sc | g |* sc ) - (tf | q |*o) = 2 i/rai. (19) 
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Finally, using (ITS]) and (|T9|) . we obtain the nuclear stiffness coefficient Cq >sc in scaling ap- 
proximation as 

_ d*AE _ m 3 
C Q , SC - -qjT - 2^2, (20) 

which differs significantly from the adiabatic one Cq^ ( fT5l) . The reasons of such a deference 
will be made clear in the next section. 



III. RESPONSE FUNCTION AND EWS FOR NUCLEAR FERMI-LIQUID 

It is necessary to make some additional assumptions for the practical calculation of the 
linear response function x{ u ) an d the corresponding EWS m^. We will restrict ourselves to 
the Landau's approximation for a nuclear Fermi-liquid and use the linearized Landau- Vlasov 
equation [ijj]. In the two-component nuclear Fermi-liquid, it is necessary to consider two 
possibilities: isoscalar excitations (when protons and neutrons move in phase) and isovector 
excitations (when protons and neutrons move in antiphase). 



A. Isoscalar excitations 

For the nuclear matter in a volume V in the case of isoscalar excitations, the linearized 
rinetic Landau- Vlasov equation has the same form as that for a one-component Fermi-liquid 

^5f + v-V r 6f - V p f eq ■ V r (6U self + U ext ) = 0, (21) 

where 5f = Sf n + Sf p = Sf(r,p;t) is the variation of the nucleon distribution (6f n for 
neutrons and 5f p for protons) in a phase space, v is the nucleon velocity, f eq = f eq ^ n + 
f eqtP = f eq {r,p) is the equilibrium distribution function, SU se if = SU se if(f,p;t) is a variation 
of the self-consistent mean field. The subscripts at V in (J2T1) indicate the variables of 
differentiation. The variation of the self-consistent field SU se if depends on the effective 
nucleon-nucleon interaction Vi n t- In the case of homogeneous nuclear matter it is given by 

SU self = I v mt (P,P) 6f{?J\ t), (22) 

where the additional factor 2 at the numerator is due to the spin degeneration. 
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The effective interaction Vi nt (p,p') is connected to the Landau's interaction amplitudes 

^ bi 

j oo 

WintCP,^) = -rpY^i^COS V)- ( 23 ) 

Here, -Pj(x) are the Legendre polynomials, # P p/ is the angle between the vectors p and and 
Np is the density of states near the Fermi surface, 

2Vp 2 df eq = Km> 
(2tt^) 3 <9e p P tt 2 ^ 3 

where e p = p 2 /2m*, m* is the effective mass of a nucleon (the definition of m* is given below), 
and pf is the Fermi momentum. In (1241) . we have used the equilibrium Fermi distribution 
function f eq = 9(ef — £ P ), where 6(x) is the Heaviside step function and Sp = p F /2m* is the 
Fermi energy. The presence of components with i ^ in sum (T2"3"j) caused by the dependence 
of the nuclear forces on the nucleon velocities. Below, we will restrict ourselves to the most 
important case where 

F ^0, Fi^O, F,> 2 = 0. (25) 
Note, that the interaction amplitude F\ determines the effective mass of a nucleon [14J] 

m* = (1 + F 1 /3)m. (26) 

We will introduce a variation of the nucleon density bp = Sp(f, t) and the isoscalar 
velocity field u = u(r,t), which are connected to a variation of the distribution function 
6f = f — feq = 8f(r,p] t) by the relations 

f 2dp _ 1 f 2dp p _ 1 f 2dp p 

where 

p = p(?,t) = J (|^s/( f 'P;*)' ^ = MO = / (|^ /e ^ f '^ (28) 

is the nucleon density. The velocity field u and the variation of the nucleon density Sp satisfy 
the continuity relation 

d 

—5p + Vpu = 0. (29) 



To check this relation, we will calculate the zero- moment of the kinetic equation (12TT) . Mul- 
tiplying Eq. (EH) by 2dp/ (2irh) 3 and integrating over p, we obtain 



—> ffi f 2df) — » — * 

di 6p + Vr ^ P 3 + J j2^f feqVr ' Vp6Uself = °- (30) 



Using Eqs. (1231) -([25]) and ([28]), we have 

V 



<)U se i f = — I F 5p+ ^-mpp- u 
N F V P F 



(31) 



Substituting Eq. (13T]) into Eq. (!30|) and taking the definition of m* ( |26|) into account, we 
derive the continuity equation ( 1291 ). 

To solve the kinetic equation (I2T1) . we assume that the external field is given by a plane 

wave 

Xoe i(q.r-ut) Then solution Q f Eq _ (J2JJ) 

can be presented as Q 

dfeg 



5f = 5Mr,p-t) 



de. 



-Vq(p) e 



i(q-f—u>t) 



(32) 



where f<f(p) is the unknown function. Substituting Eq. (1321) into Eq. (l2Tj) . we obtain the 
following equation for fq(p) 

2Vdf 



(oo - q- v) Vq\p) + q ■ v 



dfe 

l3 v int {p, p) -^T v q(P) + A o q ■ v = 0. 



(33) 



(2vr^ 

Let us expand the function v^(p, t) in a power series in the multipolarity / of a Fermi surface 
distortion 

oo 

V <£P) = ^ P l( COsd pq) V h (34) 
1=0 

where 9 pq is the angle between the vectors p and q. Using Eqs. (1231 . (I2"3|) . (I3"4"j) and (133|) . we 



15|: 



obtain the infinite set of equations for the amplitudes v\ 

+ (21 + §T^I F ^' " A o(2Z + l)Q l0 (s) = 0. 



Here, s = u/qvp and 



Qh'(s) = I I dxPi{x)-^—Pi,{x). 
2 J x — s 



With regard for condition (1251) . ( 1351) yields 

g M (s) (1 + ^/3) 



1 + F, 



where we have used the relations 15] 



n 



+ Qoo(s)(F + F F 1 /3 + F lS i) 



Qu(s) = sQ w (s) + -. 



(35) 



(36) 



(37) 



(38) 



Qio(s) = sQ o(s), 

The Legendre function of the second kind Qoo( s ) can be calculated by the use of Eq. (1361) . 
Taking the additional condition of analytical extension of Qoo(s) into the complex plane s 
into account, we can represent the function Qoo(s) as 



Qoo(s) 



l + -ln 



8 + 1 



71 



+ i-s9(l-\s\). 



(39) 




FIG. 1: Dependences of the strength function Imx(s) on the dimensionless parameter s for isoscalar 
excitations: the left panel is for the Landau damping regime, — 1 < Fq < 0; the right panel is for 
the zero-sound regime, Fq > 0. 

Let us evaluate the density-density response function assuming q = e~ l ^ r in Eqs. dS])-(j3j). 
Using the definition of the linear response function x(a>) from (jSJ) and the relations (1271) . 
(1321). (IMj) . we obtain 

= & = / */ ^ e ^«« '> = (40) 

Finally, taking (1371) into account, we obtain the density- density response function as 

XH = : ^°°ff M (41) 
1 - k(s)Q 00 {s) 



where 



1 („ , F x 2 



K ^ S) = ~W F [ Fo+ l + A/3 8 )? Qoo(s) = ^Q 00 ( s ). 



Function (HT1) has the same form as the collective linear response function in the general 
theory of collective motion (see e.g., [18]). The quantity Q 00 (s) is the intrinsic response 
function, and k(s) plays the role of the effective interaction parameter. 
In FIG. [U we present the dissipative response function 

X M H = ImxM, (42) 

which is obtained from Eq. (fldl) for two regimes: the Landau damping regime —1 < Fq < 0, 
left panel, and the zero-sound regime F > 0, right panel. Note that the zero-sound mode 
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is dumped at s < 1 (the Landau damping [16J). Here, the zero-sound wave propagates 
in phase with some particles and the energy transfer averaged over time from the wave to 
particles can be positive. The non-dumped sound wave exists in a Fermi liquid under the 
assumption s > 1 only. The dimensionless velocity of a sound wave s is determined by the 
Landau's dispersion equation [ijj 

1 - k(s)Q 00 (s) = 0. (43) 

If the dispersion equation (j4"3~l) is satisfied, both the response function (|4"T|) and the sound 
wave amplitude grow to infinity, and the sound wave propagation cannot be described within 
the framework of a linear response theory. The analysis showed jbll that the solution 
to Eq. (|43p exists (for real values of s) at F > only. This is illustrated in FIG. [TJ As 
can be seen from the right panel of FIG. [TJ the isolated root of Eq. (1431) exists at s > 1 
(zero-sound) for F > only. 

It is easy to see from the dispersing equation fj4"3]) that the velocity of the zero-sound 
increases monotonically with the interaction parameter F Q . 

It is useful to consider the solution of the dispersion equation fj43l at the asymptotic 
regime s —>■ oo (or Fq — > oo). Let us use the asymptotic expansion of the Legendre function 
of the second kind 

Qoo(,)l™ = -^-^-^--. (44) 

From relations fT4"3"l) and fT4"4"|) . we find the velocity of the zero-sound wave u = svp at 
F -> oo: 

3mm* 

Formula (145p can be compared with that for the velocity u\ of the normal sound (first-sound) 
in a classical liquid 

where K is the incompressibility coefficient. For the Fermi-liquid, the in compressibility 



%| Fo ^oo = s\ Fo ^ 00 v F = \l ^^P 2 F - (45) 



coefficient is given by 
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2 

K = 6 22- (1 + F ) w 220 MeV (46) 
2m* 



and we obtain 



ui = \l— = \ l - ± ^P 2 F . (47) 
" 9m V 3mm* 
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Taking Eqs. (|45|) - fj47|) into account, we derive 

"oliib-oo = m iIf ^oo ( 48 ) 

This result means that the velocities of the zero- and first sounds in Fermi-liquid coincide 
at a significant, Fq >> 1, repulsion between the particles. 

Using the expansion of the polarization response function X W M = RexM ®, © and 
expression (EH), one can find the EWS m_i, mi and 7713 for the Fermi-liquid (see also [71] ) as 

A 9 o A 2 ^ 4 / ja\ 

m_ 1 = --, rn^h—q, ms = R__ g . (49) 

Here, we have introduced the renormalized (due to the Fermi surface distortion effect) in- 



compressibility coefficient K' = K + 24c" 77/5, see Using relations (TlOT) and (1491) . we can 
derive the average excitation energy (the energy centroids of giant isoscalar resonances) in 
the adiabatic, Ei, and scaling E 3 , approximations: 



E ^ h i^ ^= fi vl^ (50) 

Using the dispersion relation E = huq between the excitation energy of a sound wave, E, 
and the sound velocity, u, and applying Eqs. (150!) and ff26|) . we obtain the sound velocity in 
the adiabatic, Ui, and scaling, u 3 , approximations: 



5 /(9/5 + FoK (51) 
V 6mm* V 3mm* 

By comparing Eq. ( )47l) and Eq. ( IBTl) . it can be seen that the sound velocity in the adiabatic 
approximation, u\, coincides with the first sound one, ui, and that the sound velocity in the 
scaling approach, u 3 , exceeds u\ significantly. The origin of this effect is the same as in the 
case of the nuclear stiffness coefficients Cq^ and Cq jSC , see Eqs. f|T5l) and fl20l . 

To clarify the nature of this effect, we will return to the kinetic equation (|33j) and consider 
the recurrence method of its solution. For a simplification, we neglect the external field in 
( 1331) assuming Ao = and use, instead of (1341 . the following expansion of the amplitude 
Vqip, t) into a series in the multipolarity / of a dynamic Fermi surface distortion: 

00 1 

MP) = Yl U lm(q)Ylm(p). (52) 

1=0 m=-l 
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Substituting amplitude (I521) into Eq. (1331) . using the expressions (1231) and (124")) . multiplying 
then Eq. fl33|) by the spherical function Y" z ^(p), and integrating over the angles of the unit 
vector p = p/p, we obtain the following equation for amplitudes vi m : 



m - v F q ^ G 'i \ lm 9 ' P ^ m ') ^i'm' = 0- (53) 



Here q = q/q, G t = 1 + F,/(2Z + 1), 



Im 



q-p 



l'm ,S j = C(lm,l'm') = / dttpY^p) cos9 qp Y Vm ,(p) 



V(g + l)(2f + l) 



= (_ X )m^ ^ ^(Zr00|10)(«V-m'|l,rn-m'), (54) 

where (/ 1 Z 2 m 1 m 2 |/m) are the Clebsh-Gordon coefficients. We will restrict ourselves to the 
longitudinal sound waves with u^m^o = [16J]. Taking condition ( J25l) for z/; into account, 
we obtain the following chain of recurrence equations from ( 1531) : 

SZ7 oo 7=Giz/i = 0, 

1 2 

SI/io 7=GV00 ^=^22/20 = 0, 

Vo v 15 

2 3 

15 V35 
(55) 



su l0 - i v / 47^I|(// - 100|10>| 2 !, - ^y/(2l + l)(2l + 3) |(ZZ + 100|10) | 2 u l+lfi = 0, 

Under some additional assumptions, the infinite chain of Eqs. (|55j) can be cat-off to obtain 
the analytical solution. We will consider two important cases, (i) Neglecting the Fermi 
surface distortions with multipolarity / > 2 in Eqs. (1551) . we obtain the solution 

uj = -^=v F q\/G Gx. (56) 
v3 

Consequently, the sound speed is given by 



This result coincides with that for the first sound velocity u\ of Eq. 11471) . Thus, the first 
sound regime corresponds to the excitations which preserve the spherical symmetry of the 
Fermi surface and leads to a displacement of the Fermi sphere as a whole, (ii) If we consider 
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three first equations in ([55]) and neglect the Fermi surface distortions with multipolarity 
/ > 3, then the solution to Eqs. fl55|) (now closed) gives the eigenfrequency 

u = l=v F qy/(G Q + 4/5)Gi (58) 

and the sound velocity 

The sound velocity given by Eq. ( |59l) coincides with U3, obtained in the scaling approx- 
imation (151j) . Thus, the scaling approximation for a Fermi-liquid means that all lower 
multipolarities of a Fermi surface distortion up to I = 2 are taken into account. As can be 
seen from (1491) . the model independent sum mi, as it should be, does not depend on the 
nuclear interaction (Landau's amplitudes F{). However, the last statement is not correct in 
the case of specific nuclear excitations, where the sound wave occurs due to the antiphase 
motion of the neutrons and the protons (isovector vibrations). 



B. Isovector excitations 



Below we consider the isovector excitations when protons and neutrons move in antiphase. 
In this case, we rewrite the kinetic equation ([211) for the protons and the neutrons separately: 

d - 

— Sf p + v- V r Sf p - Vp/p, e q • V r (5U PiSe if + U p>ext ) = 0, (60) 
d - 

foSfn +V' V r Sf n - V p fn : e q • V r (5U n;Se lf + U n ,ext) = 0. (61) 

We neglect the Coulomb interaction and assume N = Z. The corresponding corrections are 
not important on the description of the main characteristics of isovector giant resonances. 
Subtracting Eq. (1601) from Eq. (16T1) and introducing an isovector variation of the distribution 
function 

5f' = 6f n -5f p , (62) 
we obtain the kinetic equation for the isovector excitations 

%-5f + v ■ V r 5f - Vj eq ■ V r (SU' self + U' ext ) = 0. (63) 



13 



Here, f eq is the equilibrium distribution function which is the same for both protons and 
neutrons according to the above-made assumptions 

/^dp — 7}^ — •■ — 711 * 

( 2?r ^3 feq(r,p) = pn,eg = Pp,eq = 3^3 > ^ ' pfeq{r,P) = ~^— 6 (P ~ P^)' ( 64 ) 

The variation of the self-consistent field SU' sel j- in Eq. (1631) has the form which is similar to 
Eq. (E2D 

su' self = J ^^vumsmm, (65) 

where the effective interaction v' int (p,ff) for the isovector channel in the Landau approxima- 
tion is given by [17] (see also fl23l) ) 

00 

f 1=0 

The interaction amplitudes F[ for the isovector channel differ from the analogous ones F\ for 
the isoscalar channel in Eq. (|23|) . Thus, in contrast to the amplitude F , which determine the 
nuclear compressibility modulus (see ( H6l) ). the similar isovector amplitude Fn determines 



the coefficient of isotopic symmetry C sym in the Weizsacker mass formula |l7l. Il8| 



C sym = ^e F (1 + Fq) « 60 MeV. (67) 
Below, as it was earlier done for the isoscalar channel in Eq. (1251) . we assume that 

F o V0, F [^0, F; >2 = 0. (68) 



Solving the kinetic equation ( 1631) in the same manner as Eq. ( 1211) . we find the isovector 
response function 'x'iyj) like x(u;) from ( )4Ti) as 

*V) = g ;°'f> (69) 

where 



N F V 1 + ^i'/3 

The frequencies of isovector eigenvibrations (the poles of the response function (|69|) ) can 
be obtained from the dispersion equation 

1 - K '(s)Q 00 (s) = 0. (70) 
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The EWS m_i, m\ and (149|) for the isovector excitations take the form 

4 1 A AC 
m 1 = — — — , m, = a q , m q = ft —q . (71) 



* ^ sym a 1 1 n 

Here, we have introduced the renormalized isotopic symmetry energy m = C SJ/m + 85^/15 
and the effective mass m\ for the isovector channel, 



. m 
m = , 



1 + ki 

where Kj is the enhancement factor of the sum rule which is defined by the relation 

1+KI= i±mi. (72 ) 

l + Fi/3 V ; 

Note that, in contrast to the isoscalar sum mi (see fj49l) ). the sum m! x in (1711) is not 
model independent in sense that it depends on the effective mass m' and thereby on the 
interaction amplitudes Fi and F[. It is worth nothing that the continuity equation fl2U|) for 
the isovector excitations should be modified as well. Evaluating the zero moment from the 
kinetic equation ( 1631) . we obtain (see also (127 1) -( 1291) ) 

^V + V(l + ^)/m' = 0. (73) 

Here 

P = \{Pn + Py) = |, V = $Pn - S p p = J ^ 

-" - - _ - I I 2d P P x f ~J_ [ 2d P P X f 
U Un Up pj {2^hfm° J ~p e J {2^hfm 6J - 

Finally, the EWS ( 1TB allow one to calculate the energy centroids of isovector giant 
resonances for the adiabatic, E[, and scaling, E' 3 , approximations as 



TTl-\ __ O s „ m -P,/ / Tfln ^ ^S1 




E '^^r h f-^ q ' B -\/i^MV^ (74) 

It is useful to compare relations (1T4"1) with the corresponding expressions (1501) obtained 
for the isoscalar excitations. 

IV. FINITE NUCLEI. BOUNDARY CONDITIONS 

The above-developed approach can be directly applied to the study of the dynamic prop- 
erties of the infinite nuclear matter, where the distribution function distortion 5f has the 
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form of a plane wave in the r-space. Below, we also apply this approach to the description of 
the collective excitations in finite nuclei. For heavy nuclei, one can assume a sharp nuclear 
surface 18J. Then the variation Sf of the distribution function in the nuclear interior has 
the form of the plane wave (1521) or its projections on the states with a fixed multipolarity. 



Moreover, the equation of motion must be supplemented by the boundary conditions at the 
moving nuclear surface. 

To establish the boundary conditions, we introduce the force F which is caused by a 
sound wave and applied to a unit of the nuclear surface S, as well as the surface force 
Fs which is caused by a deformation of the nuclear surface. The general condition of the 
equilibrium for all forces applied to the free nuclear surface reads 



n ■ F\ s + n ■ F s 



0. 



(75) 



where n is a unit vector in the normal direction to the nuclear surface. Equation (J75J) 
represents the boundary condition to the dispersion equations fH3|) and (1701) . 

To evaluate the force F, we calculate the first moment to the kinetic equation (I2~T1) . 
Multiplying Eq. (j2ip by 2 dp p u / ( 2 7r K] 3 and integrating over p, we obtain the Eiler equation 



in the following form 



11 



12 



2l| 



d_ 

Peq di Uu 



where 511^ is the pressure tensor. For the isovector excitations, we obtain 12] 



5U 



5a' vfl + 5P'5, 



Vjil 



(76) 



(77) 



where 



2dp 



1 

N F 



2 
3 



^ P 2 5f (f, p, t) + ^-HpegSf/ (f, t) 

1 + Fq) e F 6p' (f, t) = C sym 6p' (f, t) 



(78) 



Let us introduce the isovector displacement field x' which is connected to the corresponding 
velocity field u' through the relation 

dx'/dt = -(1 + Kl )vf. 

Taking the continuity equation (1731) into account, we find 

Sfif = Pe„V ■ X'- 
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Finally, Eq. ([78]) yields 



SP' = CsympegV • X' ■ 



The pressure tensor 8a' in Eq. (jTTjl is given by 



where 



"'' 3m* J Vl-hY' 
3 



.... — / — 3^ - P 2 )5/V F (V^; + - ~ <* . V • f ) 



1 - 



:i+^)(i+^/3) 

3s 2 



Taking ( 1771) . ( 179i) and ( IHUl) into account, we obtain 

<5II a/ 3 = (V Q X^3 + V/3x' a ) + \ CsymPeq ~ | /if- J V • X £ 



(79) 



(80) 



(81) 



(82) 



The pressure tensor 5U UfJi determines the force F which acts from the side of the sound wave 
on a unit of the nuclear surface 

F v = n„m vil . (83) 

Using Eqs. fl82|) and flHBl . we evaluate the normal component of the force F applied to the 
nuclear surface: 



n ■ F 



S = -^T v TnSH vlt 



r=R 



r C sym p eq (l + kj) - -pip V • x' + 2fi' F r^r^XvXti 



2 \ d 

Csympeqi.1 + Kj) - - fip J div£' + 2 /4t^ * X ) 



r=Ro 

(84) 



r=R 



Let us calculate the normal component n ■ F$ of the isovector surface force Fs which 
occurs in Eq. (!75|) . To find the force Fs, we notice that a shift of protons against neutrons 
creates the additional surface energy in the case of isotopic symmetry given by [2] 

1 



8Es :Sym — -Peq r 0&sym j T dS. 



(85) 



Here, r is the mean distance between nucleons (Ro = r A 1//3 ), a sym is the isovector surface 
energy which is a parameter of theory, and r is a shift of the proton sphere against the 
neutron one. In units of ro, 

r=-(R p (t) - R n (t)) = - ((R + SR^t)) - (Ro - SR^t))) = -5R 1 (t), (86) 
r r r 



where 



5R 1 (t) = R a s (t)Y w (f) 



(87) 
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The amplitude ats(t) of isovector vibrations of the nuclear surface in Eq. (IB7I) is connected 
to the corresponding amplitude x' of the displacement field in a sound wave. To establish 
this connection, we note that, for a nucleus with the sharp edge, the displacement field in 
nuclear interior has the form (see Section 6 in |l8|) 

X * = a 1 (t)±V r (j 1 (qr)Y w (r)), (88) 

where ji(x) is the spherical Bessel function. 

Evaluating the normal component of the velocity field u' by the use of Eq. (1HH1) and 
equating it to the surface velocity d5Ri(t)/dt, we obtain 

a s (t) = -a x (t) x = qR . (89) 

According to the definition of the pressure 5Ps caused by a shift of the nuclear surface (see, 
for example, the appendix to Section 6 in [18| ), we obtain the following relation from Eqs. 
(J85D and (J86D: 

6Ps = MRxlS = 3^ ^ m6Rl - (90) 
Taking into account Eqs. (1591 and (I90l . we can evaluate the normal component (n • F s ) of 
the surface force Fs in Eq. ( 1751) . The result reads 

n ■ F s = -5P S = I p f '^\ a symai (t) Y w (f) . (91) 
6 qr {L + Kj) 

Finally, from Eqs. (1751) . (1511) . flSSj) and f|9T]) we derive the following secular equation for the 
wave number q: 



I _ 2 , 2_ , ' 

'^symPeq g pp 2 ^ F 



jl (X) + 



2 , 4 
— (tip H — 

x 3 qr (l + Ki 



j[ (x) = 0. (92) 



We point out that in the classical limit of the Steinwedel- Jensen's model at a sym — > oo, the 
boundary condition ( 1921) coincides with the similar one, j[(x) = 0, in the traditional liquid 

drop model y. 

The boundary condition (1921) allows us to find the dependence of the wave number q on 
the mass number A and to evaluate the corresponding excitation energy in finite nuclei. In 
FIG. [2], we show the dependence of the energy of isovector giant dipole resonances (IGDR) 
on the mass number A obtained by the use of the explicit solution of the dispersion equation 
(1701) and EWS (1741) . For both of them, the boundary condition ( 1921) was used. As can be seen 
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FIG. 2: Dependence of the energy of the isovector giant dipole resonances on the mass number 
obtained from the dispersion equation (|70p (solid curve 2) and from EWS (|74p (dashed lines). 
Solid curve 1 is obtained from the explicit solution of the dispersion equation ([7D|) subsidized by 
the boundary condition of the Steinwedel- Jensen's model, j[(x) = 0. For all calculations presented 
in Fig. 2, we have taken the following parameters: ro = 1.2 fm, Fi = —0.64, Fq = 0.96, F[ = 1, 
Psym = 17 MeV. The experimental data were taken from j^]. 

from FIG. [2], the lowest energy of IGDR ■^m' 1 /m'_ 1 is obtained with m'_ x and corresponds to 
the first sound regime without the Fermi surface distortions. The account of a quadrupole 
Fermi surface deformation in the sum m' 3 shifts upward the curve y7^/m[. This is due 
to the additional contribution to the nuclear stiffness coefficient caused by Fermi-surface 
distortions. Involving the higher multipolarities of the Fermi surface distortions which are 
present in the dispersion equation (1701) leads to the additional increase of the nuclear stiffness 
and the excitation energy fkvi- in FIG. [2j 

As was noted above, the dependence of the nuclear forces on the nucleon velocities (com- 
ponents with Fi and F[ in (1251) and (I66p . respectively) leads to the significant difference 
between the EWS for the isoscalar and isovector excitations. In particular, the consequence 
of this difference is the asymptotic behavior of the nuclear stiffness coefficient and the zero- 
sound velocity at an increase of the internucleon interactions F Q and Fq. In FIG. [3] and HI 
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FIG. 3: Dependences of the ratio of the first sound velocity to the zero-sound one on the interaction 
amplitude Fq for the isoscalar excitations for two values of the interaction constant F\. 

we show the dependence of the ratio of the zero-sound velocity to the first sound one on the 
interaction amplitudes Fo and Fq for the isoscalar and isovector excitations. The feature 
of the isoscalar excitations is the fact that the increase of the nucleon-nucleon interaction 
leads to a shift of the zero-sound velocity towards the first sound one (see FIG. E]). This 
means that the influence of Fermi surface distortions on the collective motion in the nuclear 
Fermi-liquid becomes negligible on the increase of the nucleon-nucleon interaction. The be- 
havior of the isovector zero-sound velocity u' is qualitatively different (see FIG. H]). With 
increase in the nucleon-nucleon interaction, the velocity u' tends to the asymptotic limit 
which significantly exceeds the corresponding first sound velocity. This is a consequence 
of the general enhancement effect of collectivity of the isoscalar zero-sound caused by the 
dependence of the nuclear forces on the nucleon velocity (see also f!72|) ). 

The enhancement factor ki for the isovector excitations denned in Eq. (172|) depends 
on the interaction constants F\ and F[. Whereas the isoscalar constant F\ related to the 
effective nucleon mass m* is well studied, the isovector constant F[ is not much studied, and 
the experimental investigation of the enhancement factor kj (|72p can help for its derivation. 

The experimental derivation of the enhancement factor in the isovector EWS m[ is con- 
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FIG. 4: The same as in FIG.O but for isovector excitations. 

nected to the investigation of the nuclear absorption cross-section a a b s {uj) of 7-quanta with 
energy hio. Let us introduce the strength function S(u, q) for the density-density response 
per unit volume V. According to Eqs. (pQ), (HOj) . and (j4"2"j) . we have 

5( W ,g) = i X (d) H- (93) 

In the case of the velocity independent forces, in accordance with Eqs. ([2]), fj4~9]) and (I7T1) . 
the strength function S(u, q) is normalized by the condition 

00 

J d{hw) hwS(uj, q) = h 2 ^ q 2 p eq . (94) 


In contrast to this, in the case of isovector excitations with the velocity dependent forces, 

the normalization condition reads 

00 

J d(M hwS(u, q ) = h 2 ^(l + F[/3) q 2 p eq . (95) 


The photoabsorption cross-section a a bs(u) is connected to the strength function by the re- 
lation 

Gabsiu) = const ■ uS(u, q). (96) 
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Here, q plays the role of the momentum which is transferred to the nucleus at the absorbtion 
of a 7-quantum. The constant in Eq. (1961) can be found from the normalization condition 
for the photoabsorption cross-section cr a b s (u)). For the isovector dipole excitations and for 



the velocity independent forces, this condition reads (Reiche-Thomas-Kuhn rule) 231 ] 



oo 

f ut- s , x 27i 2 he 2 NZ 

m 1 = d{huj) a abs (uj) = — . (97) 

/ mc A 



o 



From Eqs. (USD, flMD and flHTJ), one obtains |2J] 



, s 47r 2 e 2 NZ nl . 
Oah s \u) = —~ —uS(u,q). (98) 

Cq 2 peq A 

Since the photoabsorption occurs mainly through the giant dipole resonance, the transferred 



momentum q in (19"8"|) can be taken as q = q± = 2.08/i?o |24(, that corresponds to the classical 
boundary condition j[(x) = of the Steinwedel- Jensen's model. 

In the case of the velocity dependent forces, the normalization condition S(u, q) (1941) has 
to be replaced by condition (|95l) . and the transferred momentum q = q[ has to be calculated 
by using the boundary condition (1921 . As a result, the sum rule for a a b s {u)) takes the form 
(instead of ([97])) 

oo _ 

f , / a_ \ , s 27i 2 he 2 NZ fq[Y „ 
m[ = / d{hw)a abs {uj) = ( — ) (1 + kj). (99) 



o 



mc A \qi 



In FIG. El we demonstrate the dependence of the enhancement factor rh^/rhi on the mass 
number for a number of nuclei. We point out that the exceeding of 100% of the sum rule 
m[, which is experimentally observed for the isovector giant dipole resonances, is caused 
by the dependence of the effective nucleon-nucleon interaction on the nucleon velocity. For 
the value of the isovector amplitude F[ w 1, one can adjust (on the average) the results 
of theoretical calculations of m' l (solid line in FIG. |5]) with the experimental data. The 
nonmonotonic dependence of the experimental value of m[ on the mass number A in FIG. 
is due to the shell effects which are not taken into account within the semiclassical kinetic 
theory used in this work. 

V. CONCLUSIONS 

Using the Landau- Vlasov kinetic theory, we have studied the linear response function 
and the EWS for the isoscalar and isovector excitations in heavy nuclei and the nuclear 
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FIG. 5: Dependence of the enhancement factor m^/mi of the EWS vn!-y for isovector giant dipole 
resonances on the mass number A. The results of calculations were obtained by the use of Eqs. 
([97|) and ([99]) for two values of the isovector amplitude F[ = 1.2 (solid line) and F{ = 0.58 (dashed 
line) at the constant value of F\ = —0.64; points are the experimental data from [201 ]. 

matter. An advantage of our approach is the possibility to derive the explicit analytical 
expressions and to carry out a detailed analysis for some important nuclear characteristics. 
One of the them is the nuclear stiffness coefficient. The dynamical Fermi- surface distortion 
influences significantly the formation of the nuclear stiffness. For a slow (adiabatic) nuclear 
deformation, the nuclear stiffness coefficient is derived by the low-energy sum m_i and coin- 
cides with the stiffness coefficient of the classical (non-Fermi) liquid. This stiffness coefficient 
(adiabatic incompressibility K) causes the propagation of the first sound in a Fermi-liquid 
which is not accompanied by the Fermi surface distortions [see (fTUj) . (|47|) . and (149]) ]. In the 
general case of fast motion, the derivation of the nuclear stiffness coefficient requires to solve 
the dispersion equation (1431) . A specific role is played here by the scaling approximation. 
In this last case, the quadrupole Fermi surface distortion is only taken into consideration. 
We have shown that the stiffness coefficient in the scaling approximation is determined by 
the high-energy sum 7713 and exceeds significantly the adiabatic incompressibility K. At the 
same time, the sound velocity approaches that of the Landau's first sound [see PHI) , (loll . 
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and ([59])]. 

In the presence of the velocity dependent nuclear forces, the EWS for the isoscalar and 
isovector excitations are significantly different [see (j49p and (ITTj) ]. First of all, the EWS mi, 
which is model independent for the isoscalar excitations, becomes model dependent in the 
case of the isovector excitations. Another consequence of the mentioned difference of sums 
(f49|) and (ITTj) is the different asymptotic behavior of the zero-sound velocity on an increase 
of the nucleon-nucleon interaction (see FIG. [3]andH]). 

A feature of the isoscalar excitations is that the zero-sound velocity approaches the first 
sound one with increase in the internucleon interaction. This means that the influence of 
the Fermi surface distortion on the isoscalar collective motion in the nuclear Fermi-liquid 
becomes negligible on the increase of the internucleon interaction. In contrast to this, the 
increase of the internucleon interaction for the isovector mode leads to the asymptotic zero- 
sound velocity u' which exceeds the relevant first sound velocity. The above-mentioned 
difference between isoscalar and isovector EWS allows one to explain the fact that, in many 
cases, the experimental measurement of the EWS m' l for the isovector giant dipole resonance 
gives the more than 100 % exhaustion of the corresponding sum rule. According to Eqs. ( I7T1) 
and f!72|) . the dependence of the effective nuclear forces on the nucleon velocities generates 
the enhancement factor 1 + ki > 1, which is absent for the isoscalar excitations, in the 
sum m\ for the isovector excitations (see FIG. [5]). Note that the enhancement factor f!72|) 
depends on the isovector amplitude F[. This gives, in principe, the possibility to determine 
the interaction amplitude F[ from the fit of the EWS m[ to the experimental data. 
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